Intensity correlations and entanglement by frequency doubling 
in a dual ported resonator 
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We show that non-classical intensity correlations and quadrature entanglement can be gener- 
ated by frequency doubling in a resonator with two output ports. We predict twin-beam intensity 
correlations 6 dB below the coherent state limit, and that the product of the inference variances 
of the quadrature fluctuations gives an Einstein-Podolsky-Rosen (EPR) correlation coefficient of 
Vepr = 0.6 < f. Comparison with an entanglement source based on combining two frequency 
doublers with a beam splitter shows that the dual ported resonator provides stronger entanglement 
at lower levels of individual beam squeezing. Calculations are performed using a self-consistent 
propagation method that does not invoke a mean field approximation. Results are given for physi- 
cally realistic parameters that account for the Gaussian shape of the intracavity beams, as well as 
intracavity losses. 

PACS numbers: 42.50.Lc,42.50.Dv,03.67.Mn,42.65.Ky 
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I. INTRODUCTION 

The process of parametric down conversion has been 
used widely to generate non-classical optical fields at 
the level of single photons, as well as many photon 
fields described by continuous variables^}. At the mi- 
croscopic level non-classical correlations and entangle- 
ment arise due to the possibility of converting a single 
high frequency photon into a pair of correlated lower 
frequency photons. The down conversion process can 
be used to generate so called twin beams that have in- 
tensity correlations that are stronger than would be ob- 
tained with individual coherent states 0,0, 0>IS IE IE HI < 
Twin beams have been successfully applied to sub-shot- 
noise spectroscopy |9| and quantum nondcmolition mea- 
surements [T(il |. The generation of entangled beams in 
parametric down conversion w as p redicted by Reid and 
Drummond in the late 80's pdlll2lll3| . and is crucial for 
studies of quantum teleportation and networking [Til. HH 

mmm 

The process of second harmonic generation (SHG) 
where a fundamental field at uj\ is frequency doubled 
to create a harmonic at frequency loi = 2lo\ is comple- 
mentary to parametric down conversion, and can also be 
used for creating nonclassical light. The possibility of 
using frequency doubling may be convenient for reaching 
spectral regions that are not readily accessible by down 
conversion. It is well known that quadrature squeezing 
of both the fundamental and harmonic fields occurs in 
second harmonic generation|l9l I20L l2l| . The generation 
of multibeam correlations in second harmonic genera- 
tion is less well studied than in the case of parametric 
down conversion. Calculations have demonstrated the 
existence of correlations between the fundamental and 
harmonic fields [2^, [2^, 0, E3 including entanglement 
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FIG. 1: (color online) A dual ported singly resonant cavity 
which provides two harmonic outputs. 



between the fundamental and harmonic fields [^E [2?J and 
entanglement in type II SHG in the fundamental fields 
alone[2^, E3- The possibility of nonclassical spatial cor- 
relations in either the fundamental or harmonic fields 
alone 13(1 |3ll l32l | and of entanglement in the fundamen- 
tal field 33] has also been shown in models that include 
diffraction. 

In this work we investigate the production of two 
beams at the harmonic frequency that exhibit nonclassi- 
cal intensity correlations and quadrature entanglement. 
The device we analyze is the dual ported resonator shown 
in Fig. ^ The pump beam at ui\ is resonant in the cavity, 
while the generated harmonic exits at both end mirrors. 
As we show below the harmonic outputs exhibit strong 
quantum correlations. This can be understood naively as 
follows. While traversing the crystal to the right squeez- 
ing is generated in the fundamental and harmonic fields, 
as well as correlations between the fundamental and har- 
monic. Even though all of the harmonic leaves the cav- 
ity at the right hand mirror, the intracavity fundamental 
that generates a harmonic field on the backwards pass 
through the crystal is correlated with harmonic output 
1. The intracavity field transfers the correlations to out- 
put 2, leading to a nonzero correlation between outputs 
1 and 2. 

In order to describe this process theoretically we need 
to go beyond the usual mean field description of SHG 
squeezing^, and account for variations in the fields at 
different locations inside the cavity. To do so we com- 
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bine the linearized solutions for propagation of quan- 
tum fluctuations in traveling wave SHG0, |35| with a 
self-consistent application of the cavity boundar y co ndi- 
tions as was first done by Maeda and Kikuchi[3q. It 
was shown in Ref. that the propagation model re- 
produces the mean field results at low coupling strength, 
but predicts larger quantum noise reduction at high cou- 
pling strength. Since the propagation model relies on 
a linearized description of the quantum fluctuations its 
validity may break down in the limit of very large squeez- 
ing, where the fluctuations in the unsqueezed quadrature 
are no longer small. The accuracy of the linearized model 
was studied by comparison with numerical solutions of a 
quantum phase space model[23- It was found that the 
linearized model provides accurate predictions provided 
the normalized interaction length Cp4. l35| does not ex- 
ceed 2-3. In the results presented here, using realistic 
physical parameters and Gaussian beams, we have £ < 1. 
We are therefore confident that the linearized theory pro- 
vides an accurate description of the nonlinear resonator 
for the range of parameters considered here. 

The structure of the paper is as follows. In Sec. ITTI we 
define notation and present the solutions for the quadra- 
ture fluctuations and squeezing spectra of the harmonic 
outputs. These results are used to calculate the in- 
tensity correlations in Sec. IIIII EPR correlations and 
inseparability|37j of the outputs are demonstrated in Sec. 
IIVI We compare the EPR correlations created in the dual 
ported resonator with the alternative approach of com- 
bining two separate SHG resonators at a beamsplitter in 
Sec. and give a discussion of the results obtained in 
Sec. EH 



II. PROPAGATION MODEL OF 
INTRACAVITY SHG 

The dual ported resonator of Fig. shown in more 
detail in Fig. |21where we have separated the forward and 
backward passes to create an equivalent ring resonator 
model. Since the counterpropagating beams are strongly 
phase mismatched in the crystal we can neglect any di- 
rect interaction between them. The pump beam at fre- 
quency lo\ enters through an input coupler described by 
transmission and reflection matrices Ti , Ri (defined be- 
low) . The second harmonic outputs are allowed to escape 
through input and output mirrors. The cavity is assumed 
to be resonant at the frequency lo\ of the external pump 
beam. 

The coupled propagation of the slowly varying wave 
envelopes of the fundamental (frequency w{) and second 
harmonic (frequency uj 2 ) in the crystal is described by 
the classical equations 
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The intensities of each field are given by Ii=^-riic\£i\ 2 , 
rii is the refractive index at frequency Ui, c is the speed 
of light in vacuum, and the phase mismatch is Ak = 
2k\ — k 2 where fc.; = ujiUi/c. The coupling constants are 
given by Ki = u\d/ (cni), where d is the effective second- 
order susceptibility of the crystal. 

These equations are valid for plane waves whereas 
real experiments are typically performed with Gaussian 
beams. In order to obtain the correct value for the cou- 
pling constant when using Gaussian beams we avail our- 
selves of the results of the Boyd-Kleinman theory j^] . We 
assume that the fundamental field is a lowest order Gaus- 
sian beam with waist w (1/e 2 intensity radius) that is 
symmetrically located in the center of th e crystal and in- 
troduce scaled field amplitudes Ai = i^f e crii / ' {2hwi)£i. 
The characteristic field strength £i is chosen such that 
I Ai 1 2 gives the number of photons per second at frequency 
uji carried by the Gaussian beam. The coupled equations 
for the scaled amplitudes then take the form 
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where the complex coupling coefficient is 
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where 4>h — arg(/i), L c is the length of the crystal, and 
zri — nniUJiw 2 1 \2ttc) is the Rayleigh length of the fun- 



damental beam. The single pass conversion efficiency 
is determined by Enl = , where P\,P% are the 

powers of the two fields. The Boyd-Kleinman solution as- 
sumes that the harmonic is generated with a waist that 
gives equal Rayleigh lengths for both frequencies, and 
that the attenuation of the fundamental field due to up- 
conversion is small. As shown in Appendix ^ the frac- 
tional depletion of the intracavity field due to a single 
pass through the crystal is always small for the param- 
eters considered here. For parameters where this is not 
the case it would be necessary to use a more cumbersome 
multimode theory that accounts for the coupling of dif- 
ferent spatial modes 39] which is outside the scope of the 
present work. 

To solve for the transformation of quantum fluctua- 
tions we replace the classical fields in Eqs. © by an- 
nihilation operators Aj. The propagation equations for 
Aj are the same as the classical Eqs. (J5J. These nonlin- 
ear operator equations can be solved perturbatively by 
putting Aj(z,t) = Aj(z) + dj{z,t) where Aj{z) are the 
(classical) mean fields and &j(z, t) time and space depen- 
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FIG. 2: (color online) Propagation model of singly resonant 
SHG cavity with dual output ports. 



dent fluctuation operators which satisfy the commuta- 
tion relations [&i(z, t), cij(z', t')] — Sijd(z — z')S(t — t'), 
[a,i(z,t], dj(z' ,t')] = 0. We then linearize in the fluctu- 
ation operators, and write the resulting operator equa- 
tions as equations for classical fluctuations with the re- 
placements dj(z,t) — > aj(z,t), a,j(z, t) — » a*(z,t). The 
semiclassical theory ^ expresses the expectation values 
of symmetrically ordered quantum operators in terms of 
the classical c-numbers a,j , a* . 

The solutions to the linearized propagation equations 
are known [33.l35| . We can write the solutions in the form 
x(£,i) = N(£)x(0, t) where x = (xi,X2,yi,y2) T and 



Xj(z,t) 

%(M) 



cij(z, t) + a*(z, t), 
-i[a,j(z,t) - a*(z,t)}, 



are the amplitude and phase quadrature fluctuations. 
When ^2(0) = the transformation matrix is given by 



N = 



/Nu N 12 

N 2 i N 22 

iV 33 A^34 

V JV43 AW 
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Expressions for the matrix elements Nij(Q are given in 
Appendix^] The normalized propagation length is given 
by 



C = -j=\A 1 (0)\\ K \L e . 
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Using Eqs. © we can express £ in terms of experimen- 
tally accessible parameters as 
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At this point we note that the transformation matrix 
(JIJ is only valid when k is real. In the more general case 
of complex k propagation mixes the x and y quadra- 
tures. We will limit our study to the case of real k for 



which the analytical solutions can be expressed in terms 
of simple hyperbolic functions (the more general case in- 
volves elliptic functions). We therefore wish to have k 
real which is the case when Afc = 0. Unfortunately the 
maximum value of E^l and hence £ is obtained for|38j 
Zri = 0.176L C and AkL c = 3.26 which implies a complex 
value of k. For the analysis presented below we choose 
Afc = in order to make k real. For zero phase mis- 
match the optimum value of the beam focusing corre- 
sponds to zjii = 0.325L C , which gives a value of E^l 
that is about 40% smaller than could be obtained with 
nonzero phase mismatch. We will assume these focus- 
ing conditions in all the subsequent analysis. Numerical 
results will be given for a KNbC>3 crystal, fundamental 
wavelength of Ai = 860 nm, L c — 1 cm, d = 11 pm/V, 
m = n 2 = 2.2, and Afc = for which the optimum 
focusing is w — 21.1 /im which gives En l = 0.015 W _1 . 

A self-consistent solution for the fluctuations in the 
cavity of Fig. is found by combining the transfer ma- 
trix for the crystal propagation with the effects of mirror 
reflections and transmissions as well as intracavity losses. 
We first transform to frequency domain variables defined 
by 



aj(z,Q) — / dtaj(z,t)e 



(7a) 



a*(z,-n) = [oj-(z,-n)]* = / dta*(z,t)e lUt . (7b) 



The corresponding frequency domain quadrature fluctu- 
ations are 

Xj(z, ft) = a,j(z, O) + a,j(z, — O), (8a) 
yj (z,fl) = -i[ aj (z,n)-a*{z,-n)]. (8b) 

Note that in the frequency domain the quadrature am- 
plitudes x(z,Q),y(z,Q) are complex variables. 

We introduce the 4x4 matrices for trans- 
and reflection with nonzero diagonal el- 
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R, = diag(^/l - Tij, y/1 - T 2j , y/1 - T Xj , y/1 - T 2j ), 
where denotes the power transmittance for frequency 
i at mirror j. Residual intracavity losses due to crystal 
absorption, reflections at crystal surfaces, and mirror 
losses are lumped into effective loss beamsplitters 
indicated by £3,^4 in Fig. [21 These losses are described 
by coefficients Lij for frequency Wj at position j and 
corresponding reflection and transmission matrices 
T £j = diag(- v /l - Lij, ^1 - L 2j , y/l - L ljy y/1 - L 2j ), 

K C] = diag( v /Ll~, y/Lij, ^/L^j). 

The phase shift acquired in one cavity 
round trip is represented by the matrix D = 
diag(e in / l/cl ,e in /^ 2 ,e l °/^ 1 ,e in /^ 2 ), where the cav- 
ity free spectral range is v C i = c/(2n,iL c + 2L a ), and 
2L a is the round trip length of air in the cavity. We 
have assumed the cavity is on resonance, so that the 
phase shift is an odd function of f2. This ensures that 
D can be used with the quadrature fluctuation vector 
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x(z, CI) which contains components at ±f2. If the cavity 
were detuned it would mix the x, y quadratures and the 
round trip phase would have to be calculated separately 
for a(z, CI) and a* (z, — ft). 

Vacuum noise sources enter the cavity at mirrors 1 , 2 
and through the loss ports. We describe these by quadra- 
ture noise vectors Vj = (uij(Cl), U2j(Cl), Vij(Q), V2j(Cl)) T , 



where Uij (CI) , (CI) are frequency domain amplitude and 
phase quadrature fluctuations of frequency uji at position 

j- 

Using Eqs. 1 1181) and the above definitions the self- 
consistent solution for the intracavity fluctuations to the 
right of beamsplitter 1 is 



Xl = [I - DR 1 T £4 N(C 2 )R2T£ 3 N(Ci)]" 1 D[T 1 v 1 + R 1 T £4 N(C 2 )T 2 v 2 - R 1 T £4 N(C 2 )R 2 R£ 3 v 3 - RiR£ 4 v 4 ] 

(9) 

where I is the identity matrix. The result depends on the propagation lengths £i and £ 2 which in turn are functions of 
the pump beam power and resonator parameters. Expressions for the propagation lengths in terms of experimentally 
accessible parameters are given in Appendix [S] The vectors of output quadrature fluctuations are defined as Xj = 
(Xij, X2j, Yij ;, Y2j) T , where j labels the spatial position. The outputs can be written in terms of the self-consistent 
intracavity field as 

X! = T 2 T £3 N(Ci)x 1 - R 2 v 2 - T 2 R C3 v 3 (10a) 
X 2 = T 1 T £4 N(C 2 )R 2 T C3 N(Ci)x 1 - R lVl + T 1 T C4 N(C 2 )T 2 v 2 - T 1 T £4 N(C 2 )R 2 R £3 v 3 - TiR^v*. (10b) 



Equations (|10fl are quite general and can be used to 
describe singly or doubly resonant cavities provided the 
nonlinear propagation is phase matched and the cavity 
is on resonance at both frequencies. The general expres- 
sions for the output quadratures that result from evalu- 
ation of these equations are very cumbersome. We will 
restrict ourselves to the case of a resonant fundamental, 
and complete transmission of the harmonic at mirrors 
1,2, i.e. T 2 i = T 22 = 1. We will assume that the in- 
tracavity losses £3,^4 only affect the fundamental fields 
so L23 — L24 = 0. This last assumption is not a loss 
of generality since harmonic losses can be accounted for 
at the detectors external to the cavity. Finally, since we 
will take L13 7^ which accounts for fundamental loss 
between the two passes through the crystal, we can put 
T12 = without loss of generality. With these assump- 
tions the results for the harmonic output quadrature fluc- 
tuations can be written as 

X21 = fxiuu + /13U13 + fuu u + /21U21 + / 22 M 2 2 

(11a) 

Y 2 1 = Sll«U + #13^13 + 914VU + 021U21 + Cf 22 W 22 

(lib) 

X22 = hnun + hi 3 ui3 + huuu + ^21^21 + /i 22 m 22 

(11c) 

Y 2 2 = juVll + jl3«13 + jl±Vu + J21V21 + j 22 W 22 - 

(lid) 

Explicit expressions for the coefficients /, g, h, j are given 
in Appendix [5] Except when needed for clarity we will 
in what follows suppress the dependence on z and CI for 



brevity. 

We can use the solutions (|llfl to calculate the normal- 
ized squeezing spectra of the harmonic fields at output 
port j defined by 

S XJ = \_ ,„L?,J , (12a) 



(uaiCfiJuS^n)} 
(Y 2j (Sl)Y*(n)) 

(V2l(«)«Sl( fi )>" 



(12b) 



The spectra Sxj , Syj are normalized by the input noise 
so that S < 1 corresponds to a squeezed quadrature. To 
evaluate the spectra we make the usual assumption that 
the input noise fields at different locations and frequen- 
cies are uncorrelated so that 

< Uij{Q)u* kl {fl') > = < vyityvhin') >= 8i k 8ji5{Cl - Q!). 

(13) 

Using Eq. (|13fl we can write the output squeezing spectra 
of the harmonic fields as 



Sxi = |/ll| 2 + |/l 3 | 2 + |/l4| 2 + |/ 2 l| 2 + |/ 22 | 5 



Syi = I<?ii| 2 + l5i3| 2 + |.gi 4 | 2 + |32i| 2 + |ff 22 | 



S X2 = |/lll| 2 + |^3| 2 + |/»14| 2 + |^ 2 l| 2 + |ft 2 2| 2 



S Y 2 = |jll| 2 + |jl3| 2 + b'l4| 2 + |j 2 l| 2 + |j 22 | 5 



(14a) 
(14b) 
(14c) 
(14d) 
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FIG. 3: (color online) Amplitude squeezing (left) and phase antisqueezing (center) at fi = for the parameters given in the 
text. The right hand plot shows the uncertainty product which is almost identical for the two output ports. 



It can be readily verified using the expressions given in 
Appendix [5] that when Ci — C2 — which corresponds 
to a purely linear resonator all the squeezing spectra are 
identically unity. 

Representative results for the squeezing spectra are 
given in Fig. [3] The results are shown as a function 
of the pump beam power external to the cavity for pa- 
rameters E NL i = E NL2 — .015 W -1 , Tn = .01, and 
L\z = L14 = 0.005. The dependence of C11C2 on pump 
power is shown in Appendix^ The smaller squeezing ef- 
fect on the second output beam can be attributed to the 
fact that the mean field value of the intracavity power is 
reduced due to the SHG process after the forward pass 
through the crystal. Since, as is shown in Appendix lAl 
the normalized propagation lengths are less than 0.25 at 
the highest pump power used, we can expect the lin- 
earized semiclassical results to be accurate 251. 



III. NONCLASSICAL INTENSITY 
CORRELATIONS 



We have normalized the fields such that for a coherent 
state Sf h (Q) = (\X coh (n)\ 2 ) / (\u(n)\ 2 ) = 1, thus the 
normalized variance of the detected signal due to a co- 
herent state with intensity I is just 21. 

The corresponding formula for the intensity difference 
fluctuations is 



SI- = V2I1X21 - g\/^hx. 



22- 



The variance of the intensity difference fluctuations is 
thus 



(A|«-|) s 



{\5I-\ 2 )~W- 
2I 1 {\X 2l (tt)\ 2 ) 



■ g 2 2I 2 (\X 22 ((l)\'< 



-2gy/hhC x , 



where we have introduced the correlation coefficient 
C x = (X 21 (n)X 22 *(n) + X 21 *{Q)X 22 (n)). Normalizing 
by the sum of the shot noise variance for coherent state 
outputs with the same total intensity (2Ii + g 2 2I 2 )5(0) 
we obtain 



Given the squeezing spectra we can evaluate the quan- 
tum correlation of the intensity difference or sum of the 
two output beams. Using two detectors we measure the 
intensity difference of the harmonic outputs as 

l- = h- gh 

where g is an electronic gain factor. As the intensities 
of the two harmonic outputs are not equal due to differ- 
ent propagation lengths (1 , Qi as well as the possibility of 
unequal detector sensitivities we introduce an electronic 
gain parameter g that can be adjusted to minimize the 
noise of the intensity difference or sum. We can express 
the variance of the detected intensity difference fluctua- 
tions in terms of the squeezing spectra calculated above 
as follows. The fluctuations of a beam with intensity / 
are <5/(f2) = \/2IX(Q). The variance of the fluctuations 
is given by 



(A\SI\) 2 



(\SI\ 2 )-\(SI)\ 2 

2i(\x(n)\ 2 ). 



(15) 
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norm 



hSxi + g 2 hSx2 
h + g 2 h 



y/hhCx /, „x 
9 -j— 2T (16) 
h + g 2 h 



where Cx — Cx/S(0). When the second term is negative 
the variance can be less than unity which represents a 
nonclassical twin beam correlation. 

The optimum value of g which minimizes the fluctu- 
ations is found by putting d(A\5I-\) 2 lorm /dg = which 
gives 



9opt — 



Sx2 ± \J (Sxi — Sx2) 2 



C x 



As shown in Appendix 1X1 the intensity ratio is given by 
Ii/I 2 = ei/e 2 which can be determined from Eqs. (|Alt - 
IA2|) . For input powers up to a few Watts I\ ~ I 2 and 
S x 1 ~ S x2 . Thus the optimum g values are <? op t — ±1. 

The case of g ~ 1 corresponds to an intensity differ- 
ence which has noise greater than the shot noise limit, 
while g ~ — 1 which corresponds to the sum of intensities 
gives a strong nonclassical correlation. There is a simple 
physical explanation of this effect. A positive amplitude 
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FIG. 4: (color online) Normalized fluctuations of the differ- 
ence and sum intensities at Q — using g = ±1. The inset 
shows <? pt for the intensity sum. The parameters used were 
the same as in Fig. [3] 
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FIG. 5: (color online) Normalized inferred variance of the 
harmonic outputs at Q — 0. The parameters used were the 
same as in Fig. [3] 



fluctuation in the harmonic output at port 1 corresponds 
to an increased depletion of the fundamental beam. The 
weakened fundamental then results in a smaller ampli- 
tude output of the harmonic at the second port. Thus 
the correlation function Cx is negative and the sum of 
the output intensities has a reduced expectation value. 

Using the optimum values g op t the normalized variance 
given by Eq. can be written as 



We define the normalized inference variances as 



(A\8I. 



Sxi + Sx2 ± v (Sxi 



SX2) 2 



(17) 



where the minus sign corresponds to the case g ~ — 1. 
The normalized intensity noise is shown in Fig. 01 as a 
function of pump power. The fluctuations of the intensity 
sum for g — g opt are indistinguishable from the case g = 
— 1 for the range of pump powers shown. We see that 
the nonclassical intensity correlation is stronger than the 
squeezing of each output beam shown in Fig. [3J 



IV. EPR CORRELATIONS AND 
ENTANGLEMENT 

The presence of nonclassical twin beam correlations 
motivates evaluating the presence of EPR correlations 
and entanglement between the two output beams. Op- 
tical beams with quadrature fluctuations that embody 
EPR correlations were first demonstrated in 1992|4l). 
The essence of the EPR paradox is the ability to infer 
an observable of one system from a measurement per- 
formed on a second system spatially separated from the 
first. Hence a conditional variance is used to quantify 
the degree of EPR correlation. As shown by Reid^3 a 
linear estimate of the inference variance can be used as a 
sufficient condition for the presence of the EPR paradox. 



(AX)? nf 
(AF)f nf 



(1*21 

(\Y 21 - 



-gxX 22 \ 2 )/S(0) 
g Y Y 22 \ 2 )/S(0). 



(18) 



Here gx,gy are real gain parameters that are chosen to 
minimize the inference variances. The condition for EPR 
correlations is 



Vepr 



(AX)l { (AY)f ni 



< 1. 



The variances are individually minimized with the 
choices 



gx, op t 
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where Cy = (Y 21 (n)Y 22 *(n) + Y 21 *(n)Y 22 {n))/S(0). The 
minimum of the inference product is thus 



V, 



EPR 



\_SxiSx2 



1^2 " 
4°X. 



\SyiSy 2 



Sx 2 Sy2 



(20) 



We plot the inferred variance product in Fig. as a 
function of the fundamental pump power. The variance 
product is less than 1, implying the outputs are EPR 
correlated, for pump power above about 30 mW. 

The presence of EPR correlations are a sufficient but 
not a necessary condition for entanglement of the two 
output beams. A necessary and sufficient condition for 
entanglement of Gaussian states is the inseparability of 
the density matrix describing the two output modes. 
This can be verified using the criterion of Duan et al.|37j|. 
In our notation this criterion can be written as 



A(|a|A 21 + -X 22 ) 
a 



A(\a\Y 21 - -Y 22 ) 
a 



< 2-5(0) (« 2 + 



7 



where a is a real parameter. We have included the factor 
of 25(0) on the right hand side to account for our nor- 
malization of the commutators \Xg (t), Yki(t)] = 2iSikSji, 
which is different than that used in Ref. j33|. For the 
dual ported cavity the minimum is obtained for a ~ 1 so 
that the inseparability criterion can be written as 



Vdgcz — j (Sxi 



?X2 + Cx + Syi + Syi ~ Cy) < 1- 

(21) 



In order to facilitate comparison with Vepr defined in 
Eq. (|20|l we have included a factor of 1/4 in the defini- 
tion of Vdgcz ■ Thus entanglement is indicated for both 
criteria by V < 1. Figure [S] shows that Vdgcz < 1 for all 
finite values of the pump power, and that Vdgcz < Vepr. 
These results verify that the harmonic fields at the two 
ouput ports are always entangled, but strong EPR cor- 
relations are only present when the pump power exceds 
a threshold value. 



V. EPR CORRELATIONS ON A BEAM 
SPLITTER 

An alternative approach to creating EPR correlations 
is to combine two individually squeezed beams on a beam 
splitter [l4j[. as shown in Fig. The two input beams are 
labeled with subscripts 1,2 and the two output beams 
are labeled a, b. For notational convenience we drop the 
first subscript labeling the harmonic frequency. We can 
choose the phase of the incident beams such that the 
harmonic field fluctuations transform as 



a b = 



a\ — Mi 



V2 



(22a) 
(22b) 



The quadrature fluctuations of the output fields are 



x a 




+ Y 2 ) 


(23a) 


x b 




~Y 2 ) 


(23b) 


Y a 




-x 2 ) 


(23c) 


Y b 




-x 2 ). 


(23d) 



Following the same procedure as in the analysis of the 
dual ported cavity we define the normalized inference 
variances at frequency fl as 



(AX)f nf 

(AF)fnf 



(\X a (Q) - g x X b (Q)\ 2 
6(0) 

(\Y a (fl) - g Y Y b (n)\ 2 ) 
5(0) 



(24a) 
(24b) 



Pin 




FIG. 6: (color online) Method for generating EPR corre- 
lations by mixing individually squeezed beams on a 50/50 
beamsplitter. 



The optimum g factors which minimized the inferred 
variances are given by Eqs. (|19|l to be 



3x,o P t — 



(x a (n)x b *(n) + x a *(n)x b (n)) = Cxah 

25(0)S Xb 2S Xb 
{Y a (Q,)Y b *(Q,) + Y a *(Sl)Y b (Q,)) _ C Yab 



2S 



Yb 



25(0)S Yb 
Hence, the inference product is 

[SxaSxi, - \C 2 Xab ] [S Ya S Y b - \C Yab \ 



Vepr — 
Using Eqs. (E3 we find 

Sxa = Sxi, 

SYa = SVfc 

Cxah 
CYah 



SxhSYh 



XI -I- Sy2 



Yl -T OX2) 



Sxi ~ Sy2 

Syi — Sx2- 



When the input beams 1 , 2 are generated in equivalent 
resonators we have Sxi — Sx2, Syi = Sy2 and the EPR 
correlation reduces to 



Vepr 



<iS 2 xl S- 



Yl 



(Sxi + Syi) 2 



where the condition for an EPR paradox is Vepr < 1- 
As seen in Fig. |3| second harmonic generation results in 
nonideal squeezed states with SxiSyi = P > 1- We can 
then write Vepr = AS 2 xl p 2 / (S 2 X1 +p) 2 . 

The EPR correlation is shown in Fig. for several 
values of p. We see that for strong amplitude squeezing 
the EPR correlation is quite insensitive to the degree of 
excess noise in the phase quadrature. On the other hand 
for moderate levels of squeezing, less than about 4 dB, 
Vepr increases significantly with p. Referring to Fig. OHwe 
see that P in = 0.5 W gives S x = -2 dB and p = 2.2. For 
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FIG. 7: (color online) EPR (solid lines) and inseparability 
(dashed line) correlation coefficients for nonideal squeezed 
beams combined on a beamsplitter. 



Fig. [SJthe value of Vepr generated using two sources and 
a beamsplitter depends strongly on the characteristics 
of the source. The upper curve labeled linear resonator 
shows the results obtained by using one of the outputs 
from a dual ported resonators with the parameters used 
in Figs. |«Jlto[3] This is clearly suboptimal as only half of 
the generated harmonic light is being used, and no EPR 
paradox is seen for moderate pump powers in this case. 
The lower curve labeled ring resonator corresponds to 
an optimized ring resonator which has the same param- 
eters as for the dual ported resonator except -Eatl2 = 0. 
This resonator generates much stronger squeezing, and 
also a smaller Vepr than the dual ported resonator. We 
defer until the next section a discussion of the optimal 
approach to generating EPR correlations using SHG. 



VI. DISCUSSION 




0.2 0.4 0.6 0.8 1.0 

P in (W) 

FIG. 8: (color online) Correlation coefficient Vepr produced 
by mixing outputs from dual and single port resonators. The 
correlation is calculated as a function of the total input power, 
so that in the case of two linear or ring resonators, each one 
is pumped by Pm/2 as shown in Fig. H2 



these values the EPR correlation obtained by mixing two 
individually squeezed beams on a beamsplitter is Vepr = 
1.1, so there is no EPR paradox. On the other hand for 
Pin = 0.5 W the dual ported resonator gives Vepr = 0.7. 
Thus the dual ported resonator is able to generate much 
stronger EPR correlations even though the output beams 
have a lower level of squeezing than is required using a 
beamsplitter to mix two squeezed sources. 

When combining two equivalent squeezed sources on a 
beamsplitter the DGCZ criterion as given by Eq. i|21|) 
takes on the simple form Vdgcz = Sx , which is indepen- 
dent of the parameter p. Thus, as seen in Fig. [7| non- 
separable beams can always be created by beamsplitter 
mixing, even using nonideal squeezed sources. 

We can also compare the EPR correlation generated 
in the dual ported resonator directly with beamsplitter 
mixing as a function of total pump power. As seen in 



We have shown that SHG in a resonator with two out- 
put ports can be used to generate output beams that 
exhibit nonclassical intensity correlations, as well as en- 
tanglement. The analysis uses a propagation model for 
the quantum fluctuations that goes beyond the usual 
mean field approximation. We have shown results for 
the zero frequency noise fluctuations. The nonclassi- 
cal correlations and entanglement will degrade as the 
frequency considered is increased, with the frequency 
dependence following the usual Lorentizan form for a 
nonlinear cavity |2l|. The use of SHG to create entan- 
gled beams, as opposed to the more commonly employed 
method of parametric oscillation or amplification, may be 
advantageous in that it provides additional flexibility in 
the choice of spectral region for the entangled beams. It is 
also possible to combine two individually squeezed beams 
with a 50-50 beamsplitter to create entanglement |l4| . 
The question therefore arises as to which method, the 
dual ported resonator or two separate resonators, is most 
efficient at creating usable entanglement. This question 
has been considered for the case of degenerate vs. nondc- 
generate parametric oscillation in Ref. 01 where it was 
shown that a single nondcgcncratc oscillator was gener- 
ally preferred. 

Referring to Fig. |S|we see that the EPR correlations 
strength for a given available pump power obtainable 
with the dual ported resonator lies in between the re- 
sults obtained by combining the outputs of two separate 
dual ported, or single port ring resonators. The results 
presented were obtained for resonator parameters given 
in the text. Although differences occur for different pa- 
rameter choices we believe the results presented are rep- 
resentative. Given available nonlinear crystals which fix 
the value of Enl the most important parameters that 
can be varied are the internal losses and input coupler 
transmissions. Results have been given for round trip 
passive losses of L13 + L14 = 1%. Reducing this further 
gives larger nonclassical effects but seems unrealistic for 
current experiments. We have also chosen a low input 
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coupler transmission of T\\ — 1%. While this does not 
give the maximum possible harmonic conversion we have 
found that it is close to optimum for generating squeez- 
ing. 

Although Fig. [5] appears to show that, as far as EPR 
correlations are concerned, it is more efficient to com- 
bine two separate resonators, this conclusion may not be 
warranted. The single output resonator which generates 
the strongest EPR correlation in Fig. [5] does so by pro- 
ducing large amounts of squeezing. For the parameters 
used in the text and P; n /2 = 0.5 W the ring resonator 
is predicted to give about 6 dB of amplitude squeez- 
ing into a single harmonic output. On the other hand 
the largest amount of amplitude squeezing ever reported 
in a SHG experiment was, to the best of our knowl- 
edge, measured to be 2.4 dB, with an inferred squeez- 
ing of 5.2 dB|43(. It may therefore not be possible to 
achieve the level of squeezing predicted in the theoreti- 
cal model. The difficulties include parasitic effects such 
as blue light induced infrared absorption that become 
prominent when large amounts of the harmonic field are 
generated |44j. There are also more fundamental limita- 
tions present in high conversion efficiency SHG due to 
the excitation of competing parametric processes which 
have been shown to limit the level of harmonic squeezing 

that is attainable^]- 

Given these considerations the use of a dual ported 



We have introduced the conversion efficiencies e\ = 
^21 /-Pin, £2 = fW-Piiu P 21,^22 are the harmonic out- 
put powers after the first and second crystal passes, Pi n 
is the fundamental pump power external to the cavity, 
and the other parameters are defined in Sec. [HJ When 
Eml2 = L14 = T12 = Eq. IjAl \ reduces to the known 
result 46] for a single crystal ring cavity. 

In the two crystal case the conversion efficiencies are 
related by 



e 2 = 9^(1 - T i2) 2 (l - ii3) 2 ei(f - V^A) 2 

(A2) 

Using this result in IjAl \ we get a closed expression for 
ei that can be solved numerically. The normalized prop- 
agation lengths are then given by 




(A3a) 



C2 = J— V^E NL2 P~. (A3b) 
V n 2 



resonator which generates a strong EPR correlation even 
though the squeezing level and the power of each beam 
are not large may be advantageous compared to combin- 
ing separate SHG resonators on a beam splitter. Finally, 
we note that the dual ported configuration of Fig. ^ is 
also attractive in terms of experimental simplicity, com- 
pared to a two resonator plus beamsplitter arrangement. 
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APPENDIX A: EFFECTIVE INTERACTION 
STRENGTH 

In order to use Eqs. (|10|l to calculate output spectra we 
must first evaluate the interaction strength in each crys- 
tal. Because the fundamental is partially converted to 
the second harmonic in the first pass through the crystal 
the interaction in the second pass will be slightly weaker. 

The cavity geometry is shown in Fig. [3 After some 
simple algebra we find for the conversion efficiency of 
crystal 1 



(Al) 

I 

Numerical examples using the same parameters as in 
Sec. M(Enli = E NL2 = .015 W"\ m = n 2 = 2.2, 
T n = .01, T 12 = 0, Li 3 = L 14 = .005) are shown in Fig. 

The conversion efficiencies, second harmonic power, 
and propagation lengths are shown in Fig. [!|]as a function 
of the input pump power. The parameters were chosen 
to give strong nonclassical correlations, but are not op- 
timal for power conversion since ei,e 2 peak at quite low 
power, and the cavity is undercoupled at higher input 
powers. We see that the normalized propagation lengths 
Ci 7 C2 "C f so the linearized analysis used in the paper is 
reliable. In addition the fractional pump depletion due 
to a single pass of the intracavity field with power P c 
through the crystal is given by AP C /P C = V^i^VLi^n- 
For the parameters used the fractional depletion is al- 
ways less than 5%, so the Boyd-Kleinman analysis which 
is based on the assumption of an unaltered spatial form 
for the fundamental field is a good approximation. 



4Tu>/ENLlPin 

(\ = ^ — — — J. 

[2 — s/1 — T\\\/l — T 12 (2 — L13 — L14 — V ei-EjVLl-fm — V e 2-EiYL2-Pin)] 
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FIG. 9: (color online) Power conversion efficiency, second harmonic power and normalized propagation length, in the dual 
output cavity. 



APPENDIX B: COEFFICIENTS OF OUTPUT 
QUADRATURES 



The coefficients appearing in Eqs. I|ll|) are 



hi 

/l3 

fu 

hi 
1 22 



F 

e^Vl^^j - L 14 7V 2 i(Ci)iVii(C2) 
F 

e^ci yT=7WZ^jV 2 i (Ci ) 
F 

Ar 22 (Ci) - e^/^yi - T^VT - X^JVn(d)[JVu(Ci)JV22(Ci) - ^(CO^iCG)] 



F 



e «n/".i yT^7\TVl - ^14^21 (Ci )A r i2 (C2) 



F 



(Bl) 



0ii 

313 
5l4 
521 
522 



e in /^V7iTAr 4 3(Ci) 



G 



e' n /"d yi - TnvT^yi - £14^43 (Ci)^ 33 (C2) 

G 

e^ei yT^TTTyni^ (a ) 

G 

N44 (Ci) - e'^Vl - Tnyi - Li 3 yi - Li4iV 33 (C2)[iV 3 3(Ci)iV44(Ci) - ^34^1)^43 (Ci)] 



G 



'"/"ciy! _ rnV l _ Li 4 iV43(Cl)^34(C2) 

G 



(B2) 



/in 

ft-13 
ft-14 

h 2 i 
h 2 2 



e *n/„ cl yj^yi _ Ll3 N n ( Cl ) jv 21 (C 2 ) 



F 



F 

■"/^yr^TtTyi - Li 3 ymiVii(Ci)^2i(C2) 

F 



yi-Li 3 jVi 2 (Ci))iy 2 i(c 2 ) 

F 

jV2 2 (C2) - e^/^yr^rwi - ^i3%/r^T^jVii(Ci)[jVii(C2)iV22(C2) - ^fe^fe)] 

F 



(B3) 
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.711 

313 
314 
321 
322 



G 



G 



VI - Wl - £l 3 Vm^33(Cl)^4 3 (C2) 

G 

VI - L 13 N 34 ((i)N 43 (( 2 ) 
G 

N 44 (( 2 ) - e'^^Vl - TnVl - L13VI - Li 4 jV 3 3(Cl)[^33(C2)iV44(C2) - N 34 {( 2 )N 43 (( 2 )} 

G 



(B4) 



In the above expressions, we have introduced the defi- 
nitions 



F = l-e 



y/1 - T liy /1 - L 13 y/1 - Li4JVu(Cl)JV U (C2), 



G = 1 - e in /^ v/1 - Tnv 7 ! - ^Wl - ^AMCi)^^). 
The propagation matrix elements are[3H3 



1 - CtanhC Nl fQ = 

cosh £ ' cosh £ ' 



Mi(C) 

#21 (0 

iv 33 (C) 

N 43 (() = v^tanhC, N 44 (() = 1 ~ ( ta,nh ( 



-^(tanhC + Csech 2 C), N 22 (0 = sech 2 C, 
v2 

sechC, N 34 (() = - — (sinhC + CsechC), 
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